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^: 

^ . An SU(2) lattice gauge theory with two doublets of complex scalar fields is con- 

sidered. All continuous symmetries are identified and, using the nonperturbative 
methods of lattice field theory, the phase diagram is mapped out by direct numerical 
simulation. Two-doublet models contain phase transitions that separate qualitatively 

D , distinct regions of the parameter space. In some regions global symmetries are spon- 

taneously broken. For some special choices of the model parameters, the symmetry- 

(N 

J> ' breaking order parameter is calculated. The pattern of symmetry breaking is verified 

o. 

PsJ , further through observation of Goldstone bosons. 

in 

in 

O ■ I- MOTIVATION 

o' 

The Higgs mechanism hes at the heart of the standard model of electroweak interactions, 
/\ • where it is implemented efficiently through a single SU(2) doublet of scalar fields. The scalar 

^* ■ doublet's mass-squared term is chosen to be negative and analysis of small fiuctuations 

around the minimum of the classical scalar potential leads to the conclusion that the weak 
gauge bosons acquire appropriate masses. This standard model is generally viewed as the 
low-energy effective theory for something more complete, and in many extensions multiple 
scalar doublets appear. An old, but still useful, review of the Higgs mechanism where 
multiple doublets participate may be found in jl|. 

A different specific implementation of two scalar doublets is the inert doublet model of 
|2|(see also |3l|). Instead of taking both scalar doublets to have nonzero vacuum expectation 
values (vevs), the inert doublet model assumes that there is a phase in which one doublet 
has a vanishing vev while the other does not. This provides some additional motivation for 



examining the possible phases of two-doublet models in a more general way. 

In this work we use lattice field theory to study some features of the Higgs model with two 
doublets. Lattice field theory provides a nonperturbative method for studying non-Abelian 
gauge theories. An important difference of the lattice formulation from the continuum is 
the use of the gauge field link variable, which, taking values in the gauge group, allows 
calculations to be done without gauge fixing[4j. This has the immediate consequences that 
expectation values of non-gauge-invariant operators vanish and that the gauge symmetry 
can not be spontaneously brokenp, |6|. It also implies that physical states of the system are 
gauge-invariant composite objects^ These are features shared by lattice Higgs models and 
lattice QCD. 



Soon after lattice field theory was developed it was applied to Higgs models J8Nl2|. Lattice 
simulations of the one-doublet model were carried out extensively in the 1980s and earl} 
1990s with applications to the study of the phase diagram and basic particle spectrum 



phase transition 
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191, 



bounds on the scalar ("Higgs) mass[20] and the study of the electroweak finite-temperature 



-|24|. In addition, bounds on the Higgs boson mass have been o 



Dtained 



from simulations with Higgs- Yukawa theories that omit all gauge interactions 



An early observation 8 
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ll| was that the SU(2) lattice Higgs model with a single doublet 



of scalar fields in the fundamental representation should actually have only a single phase. 
There are regions in parameter space, sometimes named the confinement region and the 
Higgs region, which have a qualitatively different mass spectrum. In most of the parameter 
space these regions are separated by a phase transition. However, there is a corner of 
parameter space where the transition disappears and through which the confined and Higgs 
regions can be analytically connected. There is no (local) order parameter and no broken 
symmetry to distinguish the two regions. 

With regard to the spectrum, the low- lying states of the one-doublet SU(2) lattice Higgs 
model consist of a scalar singlet and a triplet of vector bosons. These states are massive 
in all regions of the parameter space. Note that the Goldstone bosons which emerge in 
an intermediate stage of the standard perturbative treatment of the Higgs mechanism and 
which are subsequently absorbed into the longitudinal components of the massive vector 
bosons do not appear in the nonperturbative lattice calculation. 



^ A view of the electroweak theory along these lines has been espoused by Frohlich et al. 



Dramatic qualitative changes may occur when additional scalar doublets are present in the 
theory. In particular, there are two regions of the phase diagram which are now completely 
separated by a phase transition throughout parameter space|33|. One might expect that 
these phases are distinguished by having different global symmetries, and if so, there will be 
corresponding order parameters. There may be regions of the parameter space where global 
symmetries are spontaneously broken and Goldstone bosons are present in the spectrum of 
physical states. In the present work, we study a gauge theory with two scalar doublets using 
numerical lattice simulations in which this scenario is realized. 

The lattice action is defined in Sec. [TTl and its continuous symmetries are discussed. 
Section Hm presents the numerical simulations used to determine vacuum expectation values 
that produce a map of the phase diagram of this two-doublet lattice Higgs model. Section HVl 
describes the methods used to search for spontaneous symmetry breaking. The symmetry- 
breaking order parameter is calculated and the Goldstone bosons that accompany each 
broken generator are identified. Section |V] contains a summary. 

II. LATTICE ACTION AND SYMMETRIES 

This study is based on an action for an SU(2) gauge theory with two complex scalar 
doublets where each doublet has its own global SU(2) symmetry^. On a spacetime lattice, 
the action can be written as 

S = J2 {^9^^^ + '^lt'^1' ^] + '^2['^2, U] + £i2[$l, $2]) , (1) 



where 



4 4 



(2) 



fj.= l u=l ^ ^ 

4 

-/€„ Y, ("^"(^ + t^)Ul{x)'^nix) + <^i{x)U^{x)<^nix + /i)) , (3) 

£i2[$i,<I'2] = 2Ai2$I(x)$i(x)<l>5(x)<l>2(x). (4) 
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This restricts the form of the aUowed quartic couphng terras. Terms allowed by gauge symmetry but 
which break the generic SU(2)xSU(2) global symmetry are excluded here although they may be present 
in phenomenological applications [1]. 



U^{x) is the gauge field and ^n{x) is a complex scalar doublet. Notice that the couplings 
Ai and A2 multiply more than just quartic terms, in contrast to common practice in the 
continuum. Likewise the normalization of the scalar fields $„(a;) follows conventions of 
lattice field theory rather than continuum conventions, and thus we show hopping parameters 
Kn instead of quadratic coefficients /x^. The classical relationship between the lattice and 
continuum notations may be found in [12]. All parameters and fields in the lattice action 
are dimensionless. 

The 4 degrees of freedom in a complex doublet, 

I a{x)+ib{x) \ 
$n(x) = , (5) 

\ c{x) + id[X) I 

can also be expressed as a matrix, 

(c{x) — id{x) a{x) + ih{x) \ 
(6) 
—a{x) + ih{x) c{x) + id{x) j 

which is a more convenient notation in some contexts. In this notation, the scalar terms in 
the Lagrangian become 

1 /I \ 2 

C-nl^nM = -Tr [(pl{x)ipnix)] + A„(^-Tr [ipl{x)ipn{x)] - ij 

4 

-K„ Y^ Tr [^l{x)Ul{x)<^n{x + /i)] , (7) 

>Ci2[(/?i,¥52] = ^Tr ^(/?l(x)y?i(x) Tr Lpl{x)ip2ix) . (8) 

This action has one local continuous symmetry, namely the gauge symmetry defined by 

U^{x) -> Rg{x)U^{x)Rl{x + fi) , (9) 

$i(x) ^ i?,(x)$i(x), (10) 

$2(x) -> i?,(x)$2(x), (11) 

where Rg{x) is an element of SU(2). The action has two global continuous symmetries, 
one for each scalar doublet, which will be referred to as intradoublet symmetries. They are 
defined by 

ipi{x) -^ ipi{x)Ri, (12) 

ip2ix) -> ip2ix)R2, (13) 



where i?i and i?2 are elements of SU(2). Finally, the action acquires an additional global 
continuous symmetry in the special case of {hi=k,2, '^i='^2='^12)- This additional symmetry 
will be called the interdoublet symmetry, and it is defined by 

*'W) ^ «J*'WV (14) 

where R12 is an element of U(2). 

It is important to understand the intricate connection between the interdoublet and 
intradoublet symmetries. To elucidate this connection, use the explicit form 

Rn=\ .^ . (15) 



-e *^"sin7„ e*""cos7„ 



for ra = 1 or 2 which gives 




<l>n{x) -^ cos7„e^""$„(x)+sin7„e^^"<l>e„(x), (16) 

<l>cn{x) -^ cos7„e-'""$c„(a;) -sin7„e-'^"$„(x), (17) 

where ^cn{x) = 'iT2$*(x). The interdoublet symmetry can be parametrized in a similar 

$i(x) \ / $i(x) 

fashion, but it is convenient to have it act on | rather than on I | , giving 

^c2{x) j \^2{x) 

gi<5i2g-iai2cQg^^2 e''^i2gi/3i2sin^^ \ /$l(x) \ 

.gi5i2g-i/3i2 sin 712 e*'^i2giai2 COS 712 j I $02(2;) j 

Now we can combine Eqs. ( !T6|) . (fT7|) and ( !T8|) to write down the transformation of our scalar 
fields under an arbitrary global transformation: 

$1 -^ cos7i2e^('^i2-"i^) (cos7ie^"^$i + sin7ie*'^i$ci) 

+ sin7i2e^(^^2+^^2) (cos72e-*"2<l>c2 - sin72e-*^2<l'2) , (19) 

$c2 -^ cos 7i2e^('^i2+"i=^) (cos 72e-*"2<l>c2- sin 726-^^2 $2) 

-sin7i2e^('^^^-''i2) (cos 716*"^ $1 + sin7ie*^^$ei) . (20) 

Finally we notice that two of the ten parameters (i.e. the a^, /3j, 7^ and 5i) are redundant. 
Let us choose the eight independent parameters to be 71, 72, 712, 

Pi = "1 + (512 - ^ , (21) 
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"12 + /3l2 



ai2 + A2 



"12 + A2 



P2 = 0:2 - <5l2 - 

^2 = /32 - ^12 - 

^ /3i2 - ai2 



When expressed in terms of the new parameters, Eqs. (IT9|) and (pOj) become 



(22) 
(23) 
(24) 
(25) 



$j ^ e*^ [cos7i2(cos7je*^^<l>j+sin7je*"'^<l>cj) +sin7i2(cos7fce ^^'=$cfc - sin7fce '"'=$fc)] , 

(26) 

$ci ^ e^'^ [cos 712 (cos7je"'''^<l>cj - sin7je"'"'^<l>j) - sin 712 (cos7fce''"=<l>fc + sin7fce^'"^$cfc)] , 

(27) 



where (j, /c) = (l,2) or (2,1). We now recognize the continuous global symmetries as 
SU(2) X SU(2) X U(l) X U(l) if (ki = K2, Ai = A2 = A12) is true. 



(28) 



where the parameters of the four factors are respectively (pi, a;i,7i), {p2i^2il2)i 7i2 and 6. 
Of course whenever {ni=H2i Ai=A2=Ai2) is not valid, the continuous global symmetries are 
just the intradoublet ones. 



SU(2) X SU(2) if (ki = K2, Ai = A2 = A12) is not true. 



(29) 



which amounts to choosing 712 = 0. The parameter 9 is then redundant, so neither of the 
U(l) symmetries remains whenever {k,i=k,2i Ai=A2=Ai2) is not valid. 

To conclude this section, return to the defining action of Eq. ([T]) and consider the special 
case of fixed-length scalar fields. 



^\{x)^l{x) = <l>^(x)<l>2(x) = 1 . 



(30) 



In this limit, the theory is independent of parameters Ai, A2 and A12. With only the Ki 



parameters remaining, the fixed-length theory 



theories and has been studied in some detail 
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Dears a notable resemblance to QCD-like 



38 



39|. 



III. THE PHASE DIAGRAM 

In numerical simulations, each scalar or gauge field is evaluated using a combination of 
heatbath and over-relaxation updates combined with an accept-reject step that accounts for 
non-Gaussian terms in the action. The algorithm contains a parameter that is tuned to 
produce a good acceptance rate. Details of the algorithm can be found in |33|; for more 



extensive discussions of algorithms see 



35 



-|37|. 



Phase transitions are readily identified, on a lattice with A^ sites, by scanning through 
parameter space and computing simple observables such as the average plaquette 

-^ J2 Trf/^(x)f/.(x + fi)Ulix + u)Ulix) , (31) 



2N 



Polyakov loops 



2N 

X n 

gauge-invariant links, where (i, j) = (1,1) or (1,2) or (2,2) 



^5^Trnf/4(x + n4), (32) 



^^i ^ ^ E {^\{^)U,{x)^A^ + /i) + h.c.) , (33) 



and the mixed vev 



^5^|<l>I(x)<l>2(x)|'. (34) 

X 

Note that the mixed gauge- invariant link {i 7^ j) and mixed vev do not preserve the intradou- 
blet symmetries, and recall that Polyakov loops are order parameters for confinement in the 
pure gauge theory and are sometimes used to provide a nonrigorous definition of confinement 
in the one-doublet SU(2)-Higgs model J40|. Figs. [Hand |2] show examples of scanning through 
Ki = K,2 values while holding /3 fixed in the fixed-length theory. For the mixed invariant link 



and mixed vev, statistical errors (not shown) scale inversely with V A^-y/^configurations at 



small K but they scale inversely with -^/^^configurations at large n. This behavior is indicative 
of the spontaneous breaking of intradoublet symmetries. The Polyakov loops are affected 
by the phase transition at /3 = 4 but not at /3 = 0.25, suggesting that the phase transition 
separates a confinement region from a Higgs region at large /3 only. This is precisely how the 
confinement /Higgs transition melts away in the one-doublet SU(2)-Higgs model as well|8|]. 
Figs. [Hand |2] also indicate that the location of the phase transition is rather insensitive to 
the size of the lattices employed. The average plaquette and gauge-invariant link undergo a 
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Figure 1: Indications of the phase transition from a variety of observables for the fixed-length theory 
at /? = 4.0. Three different lattice sizes are shown. 

qualitative change at the phase transition for all /3, and in practice the gauge-invariant link 
is a convenient first diagnostic when searching for the phase transition. 

The phase diagram for the fixed-length theory is shown for three different /3 values in 
Fig. 131 Since large /3 corresponds to weak gauge coupling, it is not surprising that one finds 
two orthogonal phase transitions: one separating the Higgs and confinement phases of the 
first scalar field (and therefore essentially independent of ^2 in the figure) and the other for 
the second scalar field (essentially independent of Ki). This divides the k,i,K2 plane into 
four regions but these are not four separate phases as is evident from results at smaller /3. 
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Figure 2: Indications of the phase transition from a variety of observables for the fixed-length theory 
at /? = 0.25. Three different lattice sizes are shown. 

At /3 = 2 the phase transitions affect one another near their mutual crossing point, and for 
/3 = 1 only a single phase transition is evident. The corresponding data for Ai = A2 = 1, 
with A12 = 0, are given in Fig. HI 

As is clear from Figs. [3] and HI there is always a single phase transition in a theory with 
degenerate scalar fields (ki = K2) but nondegenerate fields typically have more. Choosing 
Ki = 2/i;2 for definiteness, it is not clear from Fig. [3] which observables display a qualitative 
change at which phase boundaries, so this information is provided in Figs. [5] and [6l The 
Polyakov loops are zero in the Rq region but with no clear transition for small (3, while the 
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Figure 3: Phase diagram for the fixed-length theory as computed for three different /? values on 16 
lattices. R12 is the broken-symmetry region where L12 of Eq. (j33p is nonzero. All other regions are 
analytically connected at small /3 but quantitatively distinguished at larger /3 by finding large Ln 
values in region Ri, large L22 values in region R2, and small values for both Lu and L22 in region 
Rq. 

observables that mix $1 and $2 display their qualitative change at the R12 boundary. These 
results suggest that the Rq region be viewed as the confinement region, and the -R12 region 
is the phase of broken intradoublet symmetry. 

All of the simulations discussed so far have used A12 = 0, but it is interesting to explore 
nonzero values of this parameter since the lattice action has an interdoublet symmetry when 
{k,i=k,2, '^i='^2='^12)- The effect on the phase diagram due to variation of A12 is plotted in 
Fig. [7] for the case of /3 = 8.0, Ai = A2 = 1- The two phase transition lines, which were 
essentially straight and orthogonal at A12 = in Fig. HJ bend toward one another at large 
hopping parameters as A12 is increased. This pinching of the phase of broken intradoublet 
symmetry continues until that phase is reduced to a single line at A12 = 1- That line runs 
along Ki = K,2 which is precisely where the extra interdoublet symmetry is manifest in the 
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Figure 4: Phase diagram for the theory with Ai = A2 = 1 and A12 = 0, as computed for three 
different /3 values on 16^ lattices. R12 is the broken-symmetry region where L12 of Eq. (j33p is 
nonzero, region i?i has large Ln, region R2 has large L22, and region Rq has small Ln and ^22- 

lattice action. 

As A12 is increased beyond 1, the phase of broken intradoublet symmetry vanishes and a 
region of hysteresis emerges, bounded in Fig. [7] by dashed lines. Only one, not both, of the 
scalar fields is in its Higgs phase in the region between the dashed lines, meaning that the 
phenomenology of this region is similar to either the Ri region or the R2 region. Which of 
these options is realized between the dashed lines depends upon how the dynamical system 
enters the region. For example, if Ki is gradually increased to pass through that region, then 
there will be no qualitiative change in our standard suite of observables as the system enters 
the region, but there will be a qualitative change as the system exists from the region (by 
crossing the second dashed line). 
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Figure 5: Indications of the phase transition from a variety of observables for the fixed-length theory 
at /? = 4.0 with ki = 2k2- Three different lattice sizes are shown. 

IV. SPONTANEOUS SYMMETRY BREAKING 



A. Qualitative features 

The simulations discussed above found large fluctuations (for sufficiently large Hi and K2 
values) for the ensemble averages of Re($|(x)$2(a^))j Im($|(a:)$2(a^)), Re($|(x)$c2(a^)) and 
lui{^l(x)^ c2ix)) . The sum of the squares of these four quantities is observed to have small 
fluctuations for all Ki, K2 values, and is close to zero for small Ki, K2 but far from zero for 
large ki,K2. 
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Figure 6: Indications of the phase transition from a variety of observables for the fixed-length theory 
at /? = 0.25 with ki = 2k2- Three different lattice sizes are shown. 

The large fluctuations become smaller when an explicit symmetry-breaking term, such as 

5£=|Tr(^I(xV2(a:)) , (35) 

is added to the theory. Simulations can be performed for various values of rj and then 
extrapolated to t] = 0. For nonzero r], one finds ( Re($|(x)$2(a^)) ) 7^ but the other three 
ensemble averages are statistically zero, and (Re(${(x)$2(a^)) ) itself approaches zero as 



?7 -> 0. 

To determine which symmetry-group generators are broken, consider how these ensemble 
averages are affected by the general symmetry transformation of Eqs. (|26|) and ( 127|) . In 



14 



1.0 
0.8 
0.6 
0.4 
0.2 



0.0„ 



- [- 


1 1 

R,2- 
Ri - 

1 1 


-Ro 

1 1 




1 

_x,j=o.7; 
- R2 


Ay 

R, " 


-Ro 

1 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 




0.2 - R, 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 '0.0 0.2 0.4 0.6 0.8 1.0 



K, 



K, 



Figure 7: Phase diagram for the theory with Ai = A2 = 1 at /3 = 8.0, as computed for six different 
A12 values. R12 is the broken-symmetry region where L12 of Eq. (f33|) is nonzero, region Ri has 
large Ln, region R2 has large L22, and region Rq has small Ln and L22- For A12 > 1, hysteresis is 
observed in lieu of symmetry breaking. 



particular, the use of 



^im($l<i>2; 



Re($l$e2: 



nm(${<l> 



c2. 



leads to 



$[$2) -^ (cos7iCOS72e^("2-"^) + sin7ism72e*(^2~''^)) /Re('$l$2 

^\^c2) -^ e"^*^sm7i2COS7i2/$^$2 -"^I"^ 

+e-2^^ [- sin^ 712 sin 71 cos -f2e'^'''^'^'^ - cos^ 712 cos 71 sin 726-^^"^+^^) 

+ sin2 7i2COS7isin72e'("i+'^2) + cos^7i2sin7i cos-f2e-''-"^+^^^ ^Re ('$1$^ 



(36) 



(37) 



(38) 



Begin with the situation where the original global symmetry was only SU(2)xSU(2). As 
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noted in the previous section, this corresponds to 712 = 0. Therefore Eqs. (137|) and (138|) are 
simply 

$1$2) -^ (cos7icos72e^("^-"i)+sin7isin72e^('^2-''i)) /Re^'I'l^s)) , (39) 

$l$c2) ^ 0. (40) 



The maximal unbroken subgroup is obtained from the case of (a;i,/3i,7i) = {02, f32,l2), 
which identifies a residual global SU(2). Therefore 

SU(2) X SU(2) -^ SU(2) if {ki = K2, Ai = A2 = A12) is not true. (41) 

Now consider Eqs. (137|) and ( 138|) in the special case where the original global symmetry 
was SU(2)xSU(2)xU(l)xU(l). Now 712 7^ 0, but again the largest unbroken subgroup is 
obtained from the case of (q;i,/3i,7i) = {(^2, (^2,12), which gives 

<l>l<l>2) ^ (Re($I$2)) , (42) 

^{^c2) -^ e-2'^sin7i2COS7i2(<l'^<l'2-<l'l<l'i\ • (43) 



As will be discussed below, lattice simulations find ($2*^2) 7^ (•^I'l'i)- Therefore 712 = 
mT/2 for some integer n, and 6 remains as a symmetry generator in the theory. As a 
consequence, the symmetry breaking in this special case is 

SU(2)xSU(2)xU(l)xU(l) -^ SU(2)xU(l) if (ki = ^2, Ai = A2 = A12) is true. (44) 

B. Symmetry-breaking order parameter 

The calculations of the previous section can give qualitative information about the phase 
diagram, but getting a quantitative estimate of the order parameter in the broken phase 
requires a different approach. The system has to be forced to choose between different 
degenerate vacua by inserting an explicit symmetry-breaking term such as Eq. (135|) into the 
theory and then studying the limit as the coefficient, 77, approaches zero. The infinite-volume 
limit should be taken before removing the symmetry-breaking term. 

Results of numerical simulations on finite size (12^ and 20^) lattices are displayed in 
Fig. m For hopping parameters below the phase transition (i.e. k, < Kc), the symmetry- 
breaking vev extrapolates linearly to zero as rj vanishes. For hopping parameters above the 
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Figure 8: A vacuum expectation value that breaks SU(2)xSU(2) is graphed as a function of its 
Lagrangian coefficient. These data are from simulations for the fixed-length theory with /3 = 8 on 
(a) 12'' and (b) 20^ lattices. 

phase transition (k > Kc), the symmetry-breaking vev appears to extrapolate to nonzero 
values, except for a bending toward zero at small rj (visible in the 12^ simulation). This 
decrease reflects the fact that there is no true spontaneous symmetry breaking in a finite 
system and is due to modes whose Compton wavelength becomes larger than the lattice size 
at small rj. As the volume is increased this effect is restricted to a smaller region near r] = 
0. 

To deal with the volume effect, we recall that the effective field theory in finite volume 
for a scalar theory in the broken phase was developed long ago 4l|, |42| and was studied 
numerically in some detail for the one-doublet model 43|. 

Doing calculations of (^Tnp\Lp2) for different volumes (from 8^ to 20^) and different t] 



and using procedures which are verified by study of the one-doublet model, we can estimate 
the infinite-volume value of the order parameter. Figure [9] shows the value of the order 
parameter for the SU(2)xSU(2) symmetry breaking in two cases: the fixed-length theory 
and the theory with Ai = A2 = A12 = 1. Using a parametrization ( iTry9|(y92 ) oc (k — k^Y for 
K > Kc we can estimate the critical k. The values of /t^ corresponding to the lines in Fig. [9] 
are 0.316 for the fixed-length theory and 0.365 for Ai = A2 = A12 = 1- 

Recall the prediction from Eq. (jUj) of an extra broken U(l) symmetry in the theory when 
Ki = K2 and Ai = A2 = Ai2. This is verified by adding a symmetry-breaking term to the 
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Figure 9: An order parameter for SU(2)xSU(2) breaking after extrapolation to ij = 0. These data 
are from simulations with /3 = 8. 

theory, and then extrapolating its coefficient to zero. In fact, the simplest way to add an 
appropriate extra term is to run simulations with Ki ^ k,2 and extrapolate the results to 
Ki = H2. An example is provided in Fig. [TUl The transition from broken to unbroken U(l) is 
found to occur at the same critical hopping parameter as the breaking of the SU(2)xSU(2). 
At first glance, the kink in the Hi = 0.48 curve of Fig. [10] may be puzzling, but comparison 
to Fig. [7] makes the interpretation clear: the kink occurs at the phase transition crossed by 
varying K2 while holding Hi fixed. 

C. Goldstone bosons 



Spontaneous breaking of any continuous global symmetry generates a Goldstone boson 
for each broken generator. The three Goldstone bosons arising from SU(2)xSU(2)— i'SU(2) 
are found to couple readily to the operators 



2 



'■Pl^2ra 



(45) 



where Ti, T2 and t^ are the standard Pauli matrices. These operators are invariant under the 
unbroken global SU(2) and under the gauge symmetry (which is never broken). Examples of 
correlators for a range of rj are shown in Fig. [11] and the corresponding squared masses are 
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Figure 10: An order parameter for spontaneous breaking of the extra U(l). These data are from 
simulations with /? = 8 and Ai = A2 = A12 = 1- 

shown in Fig. [121 The small statistical errors provide convincing evidence that the Goldstone 
boson squared mass vanishes with a linear extrapolation of 77 — ?■ 0. For comparison, the graph 
also contains results for the operator 



2 



fm 



(46) 



That operator is not an SU(2) triplet and does not couple to the Goldstone bosons, but it 
does provide evidence of a heavy scalar particle in the theory. One might wish to name the 
Goldstone bosons it a and the extra scalar boson "cr" to follow familiar notational conventions. 
In addition to the direct method for obtaining the a correlation function, the projection 



method described in 



43l | was also used. This projection is given by 
3 



a 



proj 



^ |M|2 "" 



flf2Ta 



a=l 



where, on a lattice with A^ sites, 



Mr,, 



\M\ 



^Ei- 



N 



v\'^2Ta 



Ml + Ml + Ml 



(47) 

(48) 
(49) 
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Figure 11: Correlation functions for the operators defined in Eq. (I45p . These data are from simu- 
lations on 16^ X 20 lattices with /3 = 8 in the fixed-length theory. 



A possible operator for producing the Goldstone boson associated with the breaking of 
an extra U(l) symmetry present when {k,i=k,2, Ai=A2=Ai2) is 



in. 

2 



Avi - A^2 



(50) 



Numerical simulations using this operator produced sizable statistical fluctuations as shown 
in Fig. [131 but are consistent with a mass that vanishes as r/ — )■ 0. 



V. SUMMARY 



In this paper we studied an SU(2) Higgs model using lattice field theory methods. This 
approach provides a view of symmetry breaking which is different from the one familiar from 
the usual perturbative treatment of the standard model. The difference stems from the use 
of the gauge field link in the lattice formulation. This removes the requirement of gauge 
fixing and allows all quantities to be calculated in a gauge-invariant way. The physical states 
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Figure 12: Squared masses for the Goldstone bosons obtained from Eq. (j45p and the scalar boson 
obtained from Eq. (I46p . These data are from simulations on 16'^ x 20 lattices with /3 = 8 in the 
fixed-length theory. 



of the theory are described by gauge-invariant operators which are necessarily composite. 
The lattice Higgs model is in this sense not unlike QCD. 

The SU(2) lattice Higgs model with one scalar doublet was studied long ago. Regions of 
parameter space with seemingly different physical behavior were identified by examining the 
scalar and vector particle spectrum. These were associated with confined and Higgs "phases". 
However, it was suggested that in fact the model with only one fundamental scalar doublet 
has only one phase and no symmetries, local or global, are broken. Numerical simulations 
are in accord with this expectation. The low-lying spectrum of the theory consists of a 
massive scalar boson and a degenerate triplet of massive vector bosons. 

The addition of a second scalar doublet can lead to a richer symmetry structure than in the 
one-doublet model. For the model studied in this paper the global symmetry is generically 
SU(2)xSU(2) but is enlarged to SU(2)xSU(2)xU(l)xU(l) for particular parameter choices 
which allow for symmetry under interdoublet mixing. By examining the vacuum expectation 
values of a variety of operators, the phase diagram was mapped out. The confined and Higgs 
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Figure 13: Correlation functions for the operator obtained from Eq. (I50p . These data are from 



simulations on 16^ x 32 lattices with /3 = 2, ki = 0.48 and Ai = A2 



A 



12 



1. 



regions associated with the individual doublets could be identified. When the hopping 
parameters are sufficiently large a new phase, in which there is a strong correlation of the 
two doublet fields and the global symmetry is spontaneously broken, emerges. 

The spontaneous symmetry breaking was verified by calculation of the order parameter 
for some specific choices of the model parameters. This was done by the usual procedure of 
introducing an explicit symmetry-breaking term and studying the behavior of the system as 
the volume of the simulation was increased and the symmetry-breaking term was removed. 
The presence of Goldstone bosons in the broken phase was verified by calculation of the 
correlation functions of appropriate gauge-invariant interpolating operators. In addition to 
Goldstone bosons we also find scalar states which remain massive in all phases. As in the 
one-doublet model the gauge symmetry is unbroken. 

The focus of this work was spontaneous global symmetry breaking so the question of the 
spectrum of vector bosons was not addressed. In the nonperturbative lattice approach the 
vector bosons are composite particles and are expected to be massive in all regions of the 
phase diagram. This can be confirmed by a cursory examination of the correlation functions 
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of the vector operators. However, the quantitative determination of the mass is a difficult 
problem due to the plethora of operators that can be constructed which would require careful 
study of operator mixing and also decays due to the presence of light (pseudo-)Goldstone 
bosons in the theory. Such a study might give some information about the nature of the 
theory in different regions of the parameter space. Work on the three-dimensional Higgs 
model 4^ gives some insight into the difficulty and benefit of a spectrum calculation. 
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